Effect of hole-shape irregularities on photonic crystal waveguides 
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The effect of irregular hole shape on the spectrum and radiation losses of a photonic crystal waveguide is 
studied using Bloch-mode expansion. Deviations from a perfectly circular hole are characterized by a radius 
fluctuation amplitude and correlation angle. It is found that the parameter which determines the magnitude 
of the effect of disorder is the standard deviation of the hole areas. Hence, for a fixed amplitude of the radius 
fluctuation around the hole, those effects arc strongly dependent on the correlation angle of the irregular shape. 
This result suggests routes to potentially improve the quality of photonic crystal structures. © 2013 Optical 
Society of America 



Photonic crystals (PHCs) lie at the forefront of re- 
search in photonics [1] . One of the main obstacles to the 
practical implementation of PHCs is disorder originating 
from the fabrication process. Disorder sets a severe limi- 
tation to the quality factor of high-Q PHC cavities [2-4], 
and degrades the preformance of PHC waveguides by 
inducing extrinsic radiation losses [5-8] and light local- 
ization [9] in the vicinity of the guided band edge. For 
this reason, disorder in PHCs has been the subject of 
increasingly intense theoretical [2-8, 10, 11] and exper- 
imental [3,4,12-14] studies in the last decade. Several 
theoretical works [3-5,8,10] have considered the simplest 
possible disorder model, namely circular holes with ran- 
domly fluctuating radii and/or positions, and studied ra- 
diation losses, spectral properties, and light localization 
as a function of the fluctuation amplitude. An analy- 
sis based on more realistic assumptions, clarifying the 
role of the observed deviations from a perfectly circu- 
lar shape of the PHC holes [12,15], remains a major 
challenge. Recent experiments [12] have provided strong 
evidence that the relevant parameter, that quantifies the 
effect of disorder, is the amplitude of fluctuations in the 
hole area. On the theoretical side significant progress has 
been made [2,6,7,11], in particular by providing an esti- 
mate of the disorder-induced frequency shift of the band 
edge [11], and by showing that the correlation angle of 
the hole shape - the measure of how rapidly the radius 
fluctuates around a hole - strongly affects the loss rates 
in a PHC waveguide for a fixed amplitude of the radius 
fluctuations [6]. Yet, the interplay of angle correlation 
and radius fluctuations, and the relevance of hole-area 
fluctuations, are still poorly understood. 

In this letter we address this problem by employing 
the recently developed Bloch-mode expansion (BME) 
method [10]. BME is a non-perturbative mode expansion 
that converges to the exact solution of Maxwell equa- 
tions within the desired spectral range, while requiring a 
moderate computational effort, and that allows to imple- 
ment the hole-shape model of choice in a semi-analytical 
fashion. As a prototypical system we study a Wl-defect 
waveguide, but our conclusions about the effect of irreg- 



ular hole shape have a much more general extent. We 
carry out a systematic analysis for two limiting assump- 
tions within the hole-shape model, namely fixed radius 
fluctuation (within each single hole), or fixed hole area 
fluctuation (among different holes). The result clearly 
shows that varying the correlation angle has practically 
no effect on spectrum and radiation losses in the second 
case, while it is highly determinant within the first as- 
sumption. We conclude by clarifying the mechanism un- 
derlying this phenomenology and discuss how it relates 
to the fabrication process. 

We consider a triangular lattice of cylindrical holes of 
radius R = 0.3a in a dielectric slab of thickness d — 0.5a. 
where a is the lattice parameter. The structure presents 
a "missing" row of holes (i.e. a Wl waveguide). The 
permittivity is set to £2 ~ 12.11 in the dielectric ma- 
terial and to El = 1 outside. Disorder is modeled in 
the form of fluctuations in the hole profile, given by 
R{(t>) ^ R + dR{(f>), where (f> is the polar angle relative 
to the hole center. In particular, each hole is character- 
ized by a different disorder realization (no correlation be- 
tween hole shapes is assumed) defined by random Fourier 
expansion coefficients Cm as 
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The Bloch modes of the regular Wl waveguide are 
computed with the guided-mode expansion method [16]. 
In all computations we truncate the reciprocal space of 
guided modes to Gmax = 3(27r/a). To apply BME, we 
need the Fourier transforms of the dielectric profiles of 
both the regular structure, e(r). and the disordered one, 
£'(r). We take the z-axis to be orthogonal to the slab 
plane, in which case both profiles arc piecewise constant 
along z. The only nontrivial task is then the evaluation 
of the 2-D Fourier transforms of £(p) and £'{p) at z = 
inside the slab. The former is readily computed [10], 
while the latter we compute analytically to first order in 
6R. We denote the in-plane position of the center of each 
hole by , with C running over all holes in the structure, 
the area of the structure by S, g = |g|, and the m-th 
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Bessel function of the first kind by Jm{x). Then, 



£'(g) = e(g)+ 
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Here, we always take the {Cm} coefficients for each 
hole to be Gaussian random variables, whose distribu- 
tion is fully determined by (Cm) and (|CmP)- Note that 
{Cm^o) — Oi since it is an average over complex numbers 
with a random phase. We further assume Gaussian corre- 

lation along (/), such that {SR{(f>)SR{(l)')) ~ cr^e , 
where S is the correlation angle. This sets {J2m l^'mP) = 
(SR^) = a^, and the dependence with m of the second 
moments to 
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Fig. 1: (Color online) (a): Polar and cartesian plots of a sam- 
ple hole shape for two different correlation angles, assuming 
a fixed radius fluctuation a = 0.006a; (b): The underlying 
distributions of (|CmP); (c), (d): same, assuming a fixed area 
fluctuation {5A^> = 0.004a^ 

Within this scheme, the quantities a, S and (Cq) are 
still free parameters. To set them, we consider two dif- 
ferent models. The ffi'st model consists in assuming, for 
varying 6, a given magnitude of the fluctuation of the 
hole radius, namely by setting {5R) = and {SR'^) = 
const, which sets (Co) = and a = {SR"^). Fig. f 
(a) illustrates one realization of a hole within this model 
with a = 0.006a, computed for two different values of 6. 
The second model consists instead in assuming a given 
magnitude of the fluctuations in the hole area A, thus 
setting {6 A) = and (SA^) — const. The two conditions 
determine implicitly the parameters (Co) and cr, since 



In panel (c) of Fig. I, we show a hole realization when 
those conditions arc imposed, with {SA^) = 0.004a^, for 
the same two values of S as in panel (a). It can be seen 
that in this case the magnitude of the radius fluctuations 
((5i?^) depends substantially on 5, and while the profile 
for S = 0.06(27r) appears as a reasonable representation 
of what could be expected from a state-of-the-art PHC 
[12,15], the 5 = 0.005(27r) profile displays unrealistically 
large radius fluctuations. 
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Fig. 2: (Color online) Computed radiation loss rates, (a): 
{SR^} = const model; (b): {5A^) = const model. On the 
left, the guided bands of the regular structure (solid lines) 
and the light cone (dashed line) are displayed for reference. 

In Fig. 2 we plot the radiative loss rates and fre- 
quencies of the guided modes, computed for one realiza- 
tion of a waveguide of length 256a. Panel (a) represents 
the fixed-(5i?^) model, with the four possible combina- 
tions oi 5 = 0.06(27r), 5 = 0.005(27r), a = 0.006a and 
cr = 0.002a. As can be seen, within this model, chang- 
ing d can have an effect as dramatic as that of chang- 
ing a, with an order of magnitude difference present in 
the loss rates of the modes when going from the high 
5 = 0.06(27r) to the low S = 0.005(27r). In panel (b) we 
show a similar plot but for the fixcd-((5j4^) model, with 
the same two values of S and with (SA^) = 0.004a^ and 
(SA^) — O.OOfa^. Here, in contrast, changing 6 docs not 
have a very pronounced effect, with the radiative rates 
being only very slightly higher in the low-6 case. This 
single-realization data, together with the sample profiles 
given in Fig. I, suggest that fine features in the shape of 
the holes do not contribute as strongly to the disorder 
effects as do the smooth features. 

To test this further, we computed the density of states 
(DOS) of the modes close to the band edge, by a statisti- 
cal average over 500 disorder realizations. The computed 
DOS are plotted in Fig. 3, panels (a)-(b) for the fixcd- 
((5i?^) model and panels (c)-(d) for the tixcd-{SA^) one. 
To emphasize the role of the irregular hole shape, we also 
plot the DOS with the commonly adopted [3-5,8,10] sim- 
ple disorder model of constant radius fiuctuations (cor- 
responding to the limiting case S — > 00). Oscillations in 
the high-frequency tails originate from the finite length 
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Fig. 3: (Color online) DOS histograms close to the band edge 
(marked by a dashed-dotted line), for (a.): a — 0.006a; (b): 
a = 0.002a; (c): (SA^) = 0.004a^; (d): (SA^) = O.OOla^ 



of the simulated waveguide. All the computed histograms 
show a Lifshitz tail below the band edge. Modes in this 
region are spatially localized [10,17], and slow-light prop- 
agation is consequently hindered [8, 13, 14]. From the 
DOS it appears, once more, that the correlation angle 
5 has a strong effect on the distribution only in the 
{SR'^) = const case, where the broadening of the band 
edge increases considerably for increasing S. 

In all our simulations, the sum in Eq. (2) was restricted 
to I TO I < rumax- Wc have carefully checked that all re- 
sults were well converged starting at nimax = 5. This 
remark corroborates our statement that fine features in 
the disorder are not as relevant as the smooth ones. In- 
tuitively, this is due to the typical wavelength of the 
electromagnetic modes of the PHC, which produces a 
spatial averaging of features in the dielectric profile on a 
much smaller spatial scale. 

The main result of this work can be broken down 
in two statements. First - confirming the suggestion 
made in a recent experimental study [12] - it is in- 
deed the magnitude of the fluctuations of the hole area, 
{SA^), which determines the magnitude of the disorder 
effects. This is evident from Fig. 2 and Fig. 3, where a 
{SA"^) = const model is virtually independent of S. Sec- 
ond, when {SR'^) = const is imposed instead, the disor- 
der effects are strongly dependent on S, with a small cor- 
relation angle corresponding to much smaller loss rates 
and spectral broadening. The effect on loss rates was al- 
ready indicated by a perturbative analysis [6] . Insight is 
provided by comparing Figs. 1 (b) and (d), where the av- 
erage quantities (|C,„p) are plotted for the correspond- 
ing hole models in panels (a) and (c). The dashed-dotted 
lines in panels (b) and (d) mark the convergence range 
\m\ < nrimax = 5. In the case of fixed radius fluctuation 
depicted in panel (b), the values of (|Cmp) in this small- 
to| range vary considerably when varying 5. This is due 
to the fact that the parameters a and 5 determine re- 
spectively the integral and the width of the distribution 
{\Cm\^) as a function of to. In the case of fixed hole-area 
fluctuations shown in panel (d) instead, the two curves 
differ mainly in their width, while their values in the 



vicinity of m = are closer than in the previous case. 
Thus, for fixed area fluctuations, the Cm coefhcients rel- 
evant to the results vary less as a function of 5. 

The fabrication process of a PHC does not correspond 
to either of the two limiting models considered in this 
work. The cross section of the e-beam used for lithog- 
raphy and the subsequent etching process arc likely to 
set a typical length scale for the hole shape. For PHCs 
with nominally constant hole radii, this situation is how- 
ever expected to be closer to the fixed- ((5 i?^) assumption, 
as Figs. 1 (a) and (c) also clearly suggest. We conclude 
that, while hole fluctuations are the single relevant 
parameter to quantify disorder effects, this parameter is 
only indirectly determined by the details of the fabri- 
cation process, which instead directly set the amplitude 
of radius fluctuations and their correlation angle. If any 
control over these two parameters is possible in prac- 
tice, then the best PHC quality will arise from the finest 
roughness of the irregular hole profile. This result pro- 
vides a valuable indication for further improvement of 
the quality of PHCs in view of applications in photonics. 
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